Electronic Structure of Single- and Multiple-shell Carbon FuUerenes 
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We study the electronic states of giant single-shell and the recently discovered nested multi- 
shell carbon fuUerenes within the tight-binding approximation. We use two different approaches, 
one based on iterations and the other on symmetry, to obtain the 7r-state energy spectra of large 
fuUerene cages: C240, C540, Cgeo, Cisoo, C216O and C294o- Our iteration technique reduces the 
dimensionality of the problem by more than one order of magnitude (factors of ~ 12 and 20), 
while the symmetry-based approach reduces it by a factor of 10. We also find formulae for the 
highest occupied and lowest unoccupied molecular orbital (HOMO and LUMO) energies of CgQ.„2 
fuUerenes as a function of n, demonstrating a tendency towards metallic regime for increasing n. 
For multi-shell fuUerenes, we analytically obtain the eigenvalues of the intershell interaction. 

PACS numbers: 31.15.+q, 02. 90. -Hp 



'o ■ Introduction. — The m discovery 
c/3 ' of a new simple techniquetl for the production in bulk 
pj^ , quantities of fuUerenes has triggered intensive research on 
these new all-carbon molecules and the search for other 
novel forms of carbon. Subsequently, a new type of car- 
bon structure composed of multilayered needle- like tubes 
has been discovered by high-resolution transmission elec- 
tron microscopyo. Quite recently, giant nested shells of 
onion-like fuUerenes (also called hyperfullerenesB) haye 
'"O ■ been synthesized by intense electron-beam irradiationu. 
' Essentially they consist of a composition of concentric 
spherical fuUerene cages. The innermost cage is a Ceo 
molecule which is encapsulated tw-giant fuUerenes C240, 
(^540, • • 'i one following the otheiu'El. The interlayer spac- 
ing coincides with-that for bulk graphite (3.34 A). A re- 
^ [ cent investigationEl on the stability of these brand-new 
diverse forms of carbon shows that multi-sh Specifically, 
a very recent letteil3 shows that a stability transition from 
^sj , single to multilayer fuUerenes occurs when the number of 
7—i ' atoms exceeds ~ 6000. Moreover, the concentric carbon 
CO , onion structure containing no dangling bonds provides a 
ON , challenge to graphite (comprising flat sheets of carbon 
hexagons) for the most stable form of pure carbon. 

The structures and energies for several carbon 
"onions" have recently been carried out by using 
realistic atomic potentialsQ and molecular mechanics 
^ ' calculationsB. As yet there has been little spectroscopic 
Q information available regarding the electronic properties 
O of the multi-shell fuUerenes, besides total energy calcula- 
^ , tions. The exact diagonalization of a local orbital matrix 
■ scales as (for a A^— dimensional matrix). The or- 
' thonormalization step in plane-wave methods also scales 
as A^^. These two examples illustrate typical bottlenecks 
encountered when studying one of the central thrusts of 
condensed matter theory: to compute the energetics of 
very large systems. Different suggestions on how to min- 
imize this formidable problem have recently attracted a 
lot of attention (see, e.g., Ref. 9). Here, we explore dif- 
ferent ways to reduce this difficulty by using two very 
dissimilar approaches: one based on iterations and the 
other on symmetry considerations. The iteration method 
reduces the dimensionality of the problem, N , by more 
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than one order of magnitude (factors of either 12 or 
20 depending on the chosen initial state). For instance, it 
reduces N by about a factor of 20, from 2940 to 154, for 
the C2940 fuUerene molecule. Furthermore, it can pro- 
vide some precise analytic results for the eigenvalues. By 
using a different method based only on symmetry argu- 
ments, we can reduce iV by a factor of 10. These two 
approaches give the same results for the energy eigenval- 
ues. However, several distinct features of the first method 
that make it very convenient include the following: it re- 
quires the diagonalization of a fewer number of matrices, 
each one of them smaller in size, and every one of them 
with all their elements real. 

Here we focus on a simple model of the electronic 
states of the giant single-shell fuUerenes, and multiple- 
shell carbon "onions" with the number of shells up to 
seven. The strategy used in this work includes two 
parts. First, we treat each Cgo.,i2 fuUerene as a 60-n^- 
site nearest-neighbor tight-binding model. Each site cor- 
responds to one radially directed carbon tt orbital. Since 
the TT orbitals provide the dominant contribution to the 
conducting properties of the molecules, the a orbitals 
(being mainly responsible for the elastic properties) are 
not taken into account. In spite of the slight differences 
between bond lengths and structure curvature in these 
molecules, we have assumed a uniform hopping ampli- 
tude t for all the fuUerenes under consideration. This 
assertion is motivated by the fact that results in Cgo 
changed insignificantly by considering the two different 
hopping integrals. These icosahedrally-symmetric Cgo.„2 
molecules will be studied for n — 1,2,- ••,7 (i.e.. Ceo, 

^^240, C540, C96O, Cisoo, ^2160 and C294o)- 

In the last part of this work, we use an intermolccular 
tight-binding model to account for the weak intershell in- 
teraction between the fuUerene layers and analytically ob- 
tain its eigenvalues. It is important to stress that the sim- 
ple tight-binding approximation yields reasonable quan- 
titative psults for the energy levels of each individual 
fuUerenellElilHl, especially for those near the Fermi level. 
Furthermore, the intermolccular tight-binding model has 
also been applied to the electronic states of solid Cgo- 

Single-shell giant fuUerenes. — In the description of 
the electronic states of Cgo n^ molecules, we employ 
a tight-binding Hamiltonian (Hiickel theory as, for in- 
stance, used for fuUerenes in Refs. 10-12) written as 
H = —J2<ij>^ '^t '^3- Admittedly, this is a simphfied 
Hamiltonian: the "Ising model" for the electronic struc- 
ture. However, its understanding is a convenient step- 
ping stone to the study of more complex models. Also, 
simple models are more accessible for analytical or semi- 
analytical approaches. Much more realistic and compli- 
cated tight-binding Hamiltonians can be found, for in- 
stance, in Ref. 13 and papers cited therein. For car- 
bon Ceo, we have obtained closed-form expressions of the 
eigenvalues and eigenfunctions for both tt and a states as 
well as the Green function and local densitv-of states by 
using recursion and path-integral methodalj. Here we 
focus on the giant fuUerenes. 
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To obtain the energy spectrum, we use two entirely 
different approaches, each one with its advantages. The 
first approach we_use here is our generahzation of the 
recursion methodllj, which we call the "generalized block 
recursion method" . Because of limited space, here we 
present only a brief outline of this approach. We start 
from two initial states, |uo > and \vo >, formed as 



In |uo > (I Wo >), \j > denotes the tt orbitals on a pen- 
tagon (hexagon), and \j' > denotes those on the opposite 
(i.e., j and j' are antipodes) pentagon (hexagon). Also 
V stands for parity, with the value +1 or —1. Notice 
that these two pentagon (hexagon) rings satisfy the pre- 
requisite that a fivefold (threefold) symmetry axis passes 
through their centers. Further states are then iteratively 
generated by the general recurrence relation 



H \ fm l/m > +^m+l l/rn+1 l/m-1 >i (3) 



where m designates states which are m hops away from 
the initial state. The superscript designates various suit- 
able states. On the left-hand side of Eq. (3), at most 
three terms are nonzero for each application of H on 
\fm >■ Among them, f3 and 7 can each take at most two 
values among their several possible ones. The process 
terminates after a finite number of iterations. The repre- 
sentation of H thus becomes a "block-tridiagonal" form. 
However, the dimensionality of each block is not neces- 
sarily the same. The main operational difference between 
this approach and the standard recursion method is that 
additional indices (a, j3 and 7) are needed to fully specify 
the generated states, producing very simple values (i.e., 
±1, ±2 and ^/2) for the coefficients a's and 6's. Fur- 
thermore, the standard recursion method is extremely 
cumbersome to apply for Cgo-n^ (for n > 1), while very 
convenient for only n = 1. 

The advantage of this approach lies in that it reduces 
the task of diagonalizing the extremely large original 
Hamiltonian matrix for Cgo na to the diagonalization of 
four much smaller matrices and provides precise analytic 
results for some of the eigenvalues {e.g., ±1, ±2, —3). 
For instance, in order to obtain the whole energy spec- 
trum of C240, we only need to diagonalize two (due to 
V = ±1) 14x 14 (constructed from |mo >) and two 18x 18 
(constructed from |wo >) matrices instead of a 240 x 240 
matrix. Thus, the dimensionality has been reduced by 
(factors of 12 and 20) more than one order of magni- 
tude. In general, for Cgo.„2 the dimensionality is reduced 
from 60n^ to 3n^ -I- n (5n^ — n) starting from the ini- 
tial state |uo > (l^'o >)• The relative size reduction is 




(1) 



and 
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(3n2+n)/60n2 = l/20+l/60n ~ 1/20 ( {5n'^-n)/60n'^ = 
1/12- l/60n ~ 1/12) for |mo > (ko >)• Therefore, the 
dimensionahty of the original giant fuUerene problem can 
be easily reduced by more than one order of magnitude. 
The LUMO (HOMO) energy is solved from the matrix 
constructed from |uo > (|wo >) with negative parity. 

We also use another approach, which closely follows 
the method presented in Ref. 12. This alternative and 
quite different method is based on group theory analysis. 
For instance, for C240, six 24 x 24 matrices (four of them 
with complex elements) are diagonalized to obtain the 
energy spectrum. This method reduces the dimension- 
ality of the problem by a factor of 10. For comparison 
purposes, we also use a third different method for the 
smaller cages: direct exact diagonalization. This method 
runs into limitations for very large cages. Therefore, it is 
a desirable and useful goal to be able to first reduce the 
dimensionality of the problem by more than one order of 
magnitude, in order to study very large fuUerene cages. 

The energy levels from the first two methods for var- 
ious fuUerenes are shown in Fig. 1. In Table I, we sum- 
marize the following physically important quantities for 
the giant fuUerenes: energies of the HOMO, LUMO and 
their differences {i.e., the band gaps). As the size of the 
fuUerene grows, it can be seen that the band gap becomes 
smaller. These results demonstrate the tendency towards 
metallic regime for very large fuUerenes. The best loga- 
rithmic fit formulae of the HOMO and LUMO energies 
(in units of t) for Cgo.„2 molecules with n up to 7 are 

-BhOMO = -0.67402 • „-0-72956 (4) 

and 

i^LUMO = 0.14670 -n-i-^^^^ (5) 

By taking t equal to the typical value of 2.5 eV, these 
expressions become E'homo — —1.68505 • n~°-™^^ eV 
and Ei^uMO = 0.36675 • n^i-^^^a eV. Even though these 
fits are not very accurate for small n, we believe that 
Eqs. (4) and (5) provide good approximations for larger 
values of n {n > 8) where calculations are very difficult. 
From them, the band gap can also be readily inferred. 

Multi-shell carbon fuUerenes. — In this section, we 
study the effects of the interaction between consecutive 
shells on the electronic states of the fuUerene shells. In 
the stackiiigs of the onion-like structures, the pentagons 
are alignedQ along the icosahedral directions of the twelve 
vertices, and the hexagons are stacked in a manner sim- 
ilar to the AB (or Bernal) structure of graphite. Nev- 
ertheless, and as verified in Ref. 8, the interlayer inter- 
action is mainly between adjacent shells. Each single- 
shell fuUerene is now assumed to be in its ground (lowest 
energy) state and here we focus on its HOMO. We as- 
sume an intermolecular hopping amplitude Vn between 
the pairs of fuUerenes Cgo-n^ and Cgo.(„+i)2. A tight- 
binding approximation. Hi = —^^li^Vn{\4>n >< 
4>n+i\ + \4>n+i >< 4>n\), IS then used to model the in- 
tershell interaction, where > denotes the HOMO of 
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the Cgo n^ molecule and Ng is the total number of shells. 
The V^'s can be obtained, for example, through an in- 
termolecular resonance integrals calculation. At this mo- 
ment, such calculations cannot yet be accurately done 
because of the lack of knowledge about the exact atomic 
positions. However, the Vn^s are estimated to be of the 
order of a few hundred meV. Notice that while the multi- 
plicity of the HOMO for each fullerene is five, the nonzero 
matrix element only exits between those belonging to the 
same representation. In other words, molecular states 
> and \(pn+i > in must be in the same repre- 
sentation. Thus, for a iY^-shell carbon "onion", a small 
Ns X Ns matrix (constructed from Hj), with the m-th 
row and n-th column element —Vn-iSm,n-i — VnSm,n+i, 
is sufficient to fully investigate the interaction between 
the 5Ns HOMOs. The eigenvalues for Hj can then be 
analytically solved. The results for various Ns are pre- 
sented in Table H. This model with different hopping 
amplitudes can be readily applied to the threefold de- 
generate LUMOs. 

In summary, we use two different approaches to study 
the TT-state energy spectra of fullerene cages up to C2940- 
We also find formulae for the HOMO and LUMO energies 
of fuUerenes Cgo.„2 as a function of n. These approxima- 
tions might be useful since experimentally obtained car- 
bon "onions" can have ~ 70 shells (with 294, 000 atoms 
for the 70th shell, these systems are well beyond current 
computational capabilities). We obtain the gap energy 
as a function of n and show the tendency towards metal- 
lic regime for very large fullerenes. For multi-shell car- 
bon "onions" , we analytically obtain the eigenvalues of 
the intershell interaction. The HOMO, LUMO, and gap 
energies, as well as the whole energy spectra, of these 
molecules are relevant to the several important experi- 
mental techniques which probe the spectroscopy of single 
molecules. For instance, a scanning tunneling microscope 
{e.g., Ref. 16) can be used to probe the local spectroscopy 
of fullerenes. 
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FIG. 1. The energy eigenvalues (horizontal axis) and their 
corresponding degrees of degeneracy (vertical axis) for vari- 
ous fullerenes. The energies are in units of carbon-carbon hop- 
ping integral t. The respective fivefold degenerate HOMO and 
threefold degenerate LUMO levels are indicated. For Cgo.„2 , 
there are approximately 60-n^/4 = 15-n^ distinct energy 
levels. 

FIG. 2. The HOMO, LUMO and gap energies (in units of 
eV) versus 1/n for various fullerenes Cgo.nS. 
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TABLE I. Nl = 3n^ + n {Nh = 5n^ — n) is the dimensionality of the matrix constructed from |wo > {\vq >), which provide 
the LUMO (HOMO) energies. These are significantly smaller than the dimensionality N = 60n^ of the original Hamiltonian. 

Energies of the HOMO, LUMO and band gaps for C6o.„2 fuUerenes, with 71 = 1, 2, ■ ■ ■ , 7. In the HOMO and LUMO columns, 
the first values arc iu units of f and those iu jjarentheses are the results obtained by taking t equal to the typical value of 2.5 



c\. 


Molecule 


Nl 


Nh 


HOMO 




LUMO 


Gap (cV) 


Ceo 


4 


4 


-0.618033 (- 


-1.545082) 


0.138564 


(0.346410) 


1.891492 


C240 


14 


18 


-0.436772 (- 


-1.091930) 


0.059657 


(0.149142) 


1.241072 


Cz,40 


30 


42 


-0.324099 (- 


-0.810247) 


0.033671 


(0.084177) 


0.894424 


C96O 


52 


76 


-0.255010 (- 


0.637525) 


0.022169 


(0.055422) 


0.692947 


C1500 


80 


120 


-0.209564 (- 


0.523910) 


0.015971 


(0.039927) 


0.563837 


(^2160 


114 


174 


-0.177673 (- 


-0.444182) 


0.012192 


(0.030480) 


0.474662 


C294O 


154 


238 


-0.154140 (- 


-0.385350) 


0.009690 


(0.024225) 


0.409575 



TABLE II. The energy eigenvalues of Hi for various Ni 



Ns 


Eij 


jenvalues 


2 




±Vi 


3 


and 


±^v? + vi 


4 


-A±B 


and A±B 


5 


0, ±V<S - T 


and ±V<S + T 


Here A = 

S = 


^/{V^ + Vsr + Vi/2, B 

{V? + vi + 


= V(Vi - VzY + Vil2, 
Vi + Vi)/2, and 



T = ^{Vl + Vi - Vi - Ylf + 4 ViVil2. 
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